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Low-buckled silicene is a prototypical quantum spin Hall insulator with the topological quan¬ 
tum phase transition controlled by an out-of-plane electric field. We show that this field-induced 
electronic transition can be further tuned by an in-plane hydrostatic biaxial strain e, owing to the 
curvature-dependent spin-orbit coupling (SOC): There is a .Z 2 = 1 topological insulator phase for 
biaxial strain |£| smaller than 0.07, and the band gap can be tuned from 0.7 meV for e = -1-0.07 up 
to a fourfold 3.0 meV for e = —0.07. First-principles calculations also show that the critical field 
strength Ec can be tuned by more than 113%, with the absolute values nearly 10 times stronger 
than the theoretical predictions based on a tight-binding model. The buckling structure of the 
honeycomb lattice thus enhances the tunability of both the quantum phase transition and the SOC- 
induced band gap, which are crucial for the design of topological field-effect transistors based on 
two-dimensional materials. 


Two-dimensional (2D) quantum spin Hall (QSH) in¬ 
sulator [J-Q is characterized by an insulating bulk and 
gapless edge states at its boundaries fli- These edge 
states are topologically protected from backscattering of 
non-magnetic defects or impurities due to time-reversal 
symmetry, thus providing enticing concepts for novel 
mantum electronic devices with low energy dissipation 
[a, 0] • Quantized conductance through QSH edge states 
were original reported on HgTe/CdTe 0,S , and on 
InAs/GaSb [Sj, l9| quantum wells too. 

There is an intense drive to realize QSH insulators with 
controllable quantum phase transitions and tunable elec¬ 
tronic and spin properties Elm. The intrinsic spin- 
orbit coupling (SOC) in graphene is weak 0, [iB - H^ . 
but other 2D materials that may realize this phenomena 
include honeycomb lattices of bismuth atoms on a sili¬ 
con surface [l5|, an electric-field-induced QSH phase on 
few-layer black phosphorus [l6l| . a new structural phase 
of a transition-metal dichalcogenide [13, Hi, and sil- 
icene/germanene [I^, 11, 1^24|, among others pH - l^ . 

Mechanical strain and curvature (i.e., shape) [T^, [2^ 
modify the electronic and spin properties of buckled hon¬ 
eycomb lattices, and a deep understanding of their effects 
is emerging [3(il4^ . In a classic work, the dependence of 
the SOC strength Xsoc on curvature is established to 
arise from (1) on-site spin-flips among a— and tt— bands 
due to the intrinsic coupling among spin and s— and 
p—electronic orbitals, and (2) a subsequent hybridiza¬ 
tion of the hopping tt— electron with electrons from bands 
with cr—symmetry in the presence of curvature [l^ . This 
important phenomenon lies beyond those descriptions of 
the electronic structure of silicene and other group- IV 2D 
materials that are based on tt— electrons only [l^ . El. 
Curvature is of central importance for the discussion of 
2D-based topological field-effect transistors because it 
raises the intrinsic gap, and hence the temperature at 
which these devices could operate. On the other hand, a 


controllable and tunable topological quantum phase tran¬ 
sition w ill g reatly facilitate the device fabrication and op¬ 
eration |18l |. Working on silicene, we generalize Ezawa’s 
result /X.{Ez) = 2\lEz — AsocI for the dependence of the 
energy band gap A{Ez) on the electric field (A—field) 
E 2 , in which the SOC strength Asoc is assumed to be a 
constant, into: 


A{E,,e) = 2\l{e)E,-Xsocie)\, 


( 1 ) 


where Xsoc evolves with strain due to a curvature- 
induced hybridization among s and p electrons. 

Pristine silicene is a QSH insulator with a band gap 
A{E^ = 0,£ = 0) of about 1.5 meV (18 K) [2^, |2lf. 
This band gap is tunable by an A—field E^ perpendicu¬ 
lar to the buckled atomic layer [3,[2l[23 as an electro¬ 
static potential difference is established between the two 
Si atoms in the unit cell due to their height difference 21. 
Silicene behaves as a .Z2 = 1 non-trivial QSH insulator 
below a critical Ez strength \Ec\, and becomes a triv¬ 


ial band insulator for values of \Ez\ > \Ec\ [idjlUl. In 


addition, silicene is likely to be fabricated on substrates 
3^4^ and may be subject of in-plane strain e already. 
A fundamental understanding of effects of the strain on 
the topological quantum phase transition is important 
for the design of quantum electronic devices based on 
silicene and other 2D materials Hi Em, [23. 

Calculations were carried out using density-functional 
theory [E M\ as implemented in the Quantum 
ESPRESSO code [E with the CCA-PBEsol exchange- 
correlation functional [13. The SOC was included in 
the fully relativistic Rappe-Rabe-Kaxiras-Joannopoulos 
(RRKJ) ultrasoft pseudopotential scheme [i^ with a 
nonlinear core correction. The cutoff energy in the plane 
wave expansion is 50 Ry. A Monkhorst-Pack uniform k- 
grid of 36x36 X1 is employed. A vacuum region of 20 A is 
introduced along the out-of-plane {z) direction to elim¬ 
inate spurious interactions among periodic images. The 
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FIG. 1. (Color online) The electronic properties of silicene 
are independently tnned by (a) an in-plane biaxial strain e, 
and (b) an out-of-plane i?—field Ez. Subplots (c) and (d) are 
band dispersions nnder typical values of e or Ez, respectively. 
Insets are zoom-ins of the band dispersion near the IF—point 
(the Fermi level is set to zero; note the overlapping bands for 
e = 0.00 and 0.05 on (c)). Here, we will explore the combined 
effects of e and Ez on its electronic structures. 


i?—field Ez is induced by a saw-tooth potential along the 
0 —direction (c.f., Fig. 1(b)). Biaxial strain is applied to 
the silicene lattice, and e is defined as e = (a —an)/an (see 
Fig. 1(a)). Here a and an = 3.85 A are the 

strained and unstrained lattice constants, respectively. 


There is a 1.5 meV gap opening at the IF—point that 
is induced by the SOC for e = 0 and = 0. When in¬ 
dependently applied (c.f., Figs. 1(a) and 1(b)) strain and 
i?—field lead to distinct changes on the electronic band 
dispersions, as seen in Figs. 1(c) and 1(d), respectively. 
In-plane strain preserves inversion symmetry, leading to a 
gap opening that is more apparent for compressive strain 
(see the inset of Fig. 1(c), e = —0.05), thus suggesting 
a relation among Xsoc and e that -as a matter of fact- 
has been discussed in the context of graphene some time 
ago 13| . On the other hand, E^ lifts inversion symmetry 
thus removing the band degeneracy, as seen on the inset 
of Fig. 1(d) for E^ = 1 V/nm (dashed lines). The effects 
of F/-fields on the electronic properties of silicene have 
been thoroughly discussed in the past, so we continue 
exploring the effect of strain on Fig. [5] 


There are three effects of strain on the electronic struc¬ 
ture that are conveyed by Fig. [5Ka): (i) A renormaliza- 
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FIG. 2. (a) Typical band dispersions for various magnitudes 
of e: A semicondoctor-to-metal transition occurs when |£| > 
0.07, as highlighted by the tilted (red) arrows. Insets help 
highlight the band dispersions near the Fermi level around the 
A—point, (b) Band gap A (left axis) and the height difference 
21 due to buckling (right axis) as a function of biaxial strain 
£, for an Ez = Q display a linear relationship in silicene, that 
has been explicitly indicated. 


tion of the Fermi velocity near the AT-point [s^. (ii) 
An upward shift of the valence-band maxima (VBM) 
under compressive strain, and a downward shift of the 
conduction-band minima (CBM) under tensile strain at 
the F-point, (hi) A gap opening at the A—point that is 
especially evident under compressive strain. Effects (ii) 
and (iii) are responsible for the semiconductor-to-metal 
transition, and are the focus of the ensuing discussion. 

For sufficiently large compressive strain, the VBM 
crosses the Fermi level at the F—point (c.f.. Fig. UKa) for 
e = — 0.07; this crossing has been highlighted by a tilted 
red arrow). Similarly, the CBM at the F-point nearly 
touches the Fermi level for a tensile strain of e = -1-0.07: 
Silicene undergoes a transition to a metallic phase for 
strain beyond |£| 0.07. Although arising under dif¬ 

ferent mechanisms, these effects due to strain are remi¬ 
niscent of the topological-metal phase transition that is 
induced by an electric field in few-layer black phosphorus 

M- 

We next discuss the band dispersion seen around the 
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FIG. 3. (Color online) Topological quantum phase diagram 
of silicene with respect to in-plane biaxial strain e and out- 
of-plane i?—field Ez. The vertical axis is the band gap A. 
The critical i?—fields Ec, at which there is a phase transition 
from a topological insulator into a band insulator, have been 
indicated by vertical white arrows. The metallic state has a 
zero value of A that is shown in red, and the area marked by 
QSH represents the topological spin Hall state phase. 


iiT—point in the insets of Fig. [UJa). To do so, we first 
report the SOC-induced band gap A and the buckling 
height 21 as a function of e in Fig. [2](b): A increases 
monotonically as e decreases from e = -1-0.07 down to 
e = —0.07, meaning that compressive strain enhances 
the band gap A while tensile strain decreases the gap. 
Since the energy band gap never closes in going from 
pristine silicene - a Z 2 = 1 topological insulator - to 
strained silicene with |e| < 0.07, so these systems share 
the same topological classification according to the adia¬ 
batic continuity argument for transformations of the elec¬ 
tronic bands. 

The above phenomenon can be understood from the 
change of the buckling height 21 among two silicon atoms 
due to an in-plane strain {I = l{e)) shown in Fig. [5](b). 
Buckling (and hence 1) increases with compressive strain, 
thus enhancing the overlap amon g a and tt orbitals, re¬ 


sulting in a strain-dependent SOC 


13, 


I . The band gap 
5.25Z. A positive e 


A increase is almost linear on b A 
reduces I, making the structure more planar (graphene¬ 
like), so the overlap among cr and tt orbitals decreases, 
bringing Xsoc down into its “intrinsic value” (i.e., its 
value under a zero curvature) [l^ . . However, when 

e = -1-0.07, the CBM at the F—point nearly touches the 
Fermi level (tilted arrow on subplot in Fig. 2(a)), so that 
a further increase of e induces a transition into a metallic 
phase. 

The strain-dependent SOC may be more qualitatively 


explained with the following matrix form at atom A 2lj : 
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Xsoc 
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( 2 ) 


on the basis set {|s t), \Px t), \Py t), \Pz t), |s i), \Px 4-), 
\Pv 4-)! \Pz i)}- This SOC matrix produces spin flips 
among \pz) (tt) and \px), \Py) orbitals (that belong on a— 
bands prior to the orbital hybridization) and Xsoc could 
be obtained as a fitting parameter. The SOC matrix at 
atom B Hb,soc an identical form. This intrinsic 
SOC is further enhanced by hopping [13 1. 

After indicating the dependence among A and e seen 
on the first-principles data (Figs. 1 and 2), we study the 
combined effects of strain and £1—field on the topolog¬ 
ical phase transition, and display A as a function of 
the Ez and e. A{Ez,e) has a characteristic VF—shape 
IS [30, El m as seen in Fig. 3. For a given value of 
e, A decreases to zero as \Ez \ increases. For a given mag¬ 
nitude of e smaller than 7%, the value of Ez for which 
A = 0 is known as the critical field \Ec{e)\: As the field 
increases further beyond |i?c(e)|, the energy band gap A 
reopens and increases again, but the electronic state is a 
trivial insulator because, according to the bulk-boundary 
correspondence principle, the topological phase transi¬ 
tion occurs at the critical field strength |£lc(£)| when the 
band gap closes. Fig. [3] clearly shows that the critical field 
strength \Ec{e) \ increases with compressive strain: Strain 
tunes the quantum phase transition due to a buckling- 
dependent SOC. 

Fig. [3] can also be seen as a phase diagram that indi¬ 
cates the different quantum phases that can be reached 
by jointly tuning Ez and e: Between —0.07 < e <-1-0.07, 
silicene is a QSH insulator for any applied A—field in 
the range of \E\ < \Ec{e)\. On the other hand, silicene 
changes to an ordinary band insulator with \E\ > |£lc(e)|. 
This phase diagram is clearly visible in Fig. [3] 

Finally, we discuss the critical field strength E^. In 
Fig. m the first-principles data E^ are shown as a func¬ 
tion of the theoretical predictions from Eq. (1): E^^ = 
Xsoc/I- The DFT results are nearly 10 times larger than 
the theoretical values, indicating a strong screening in sil¬ 
icene. In Ref. [23|, the field-induced band gap is found 


to be suppressed by a factor of about eight due to the 
high polarizability of silicene. Here, we show that the 
screening also significantly enhances the critical A—field 
strength required to induce a quantum phase transition. 
More precisely, a linear fit in Fig. 4 yields E^^"^ ~ 9.6El^. 

In summary, using first-principles methods, we show 
that the biaxial strain e can be utilized to tune the spin- 
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FIG. 4. (Color online) The critical S—field strength 
calculated from DFT is shown as a fnnction of the theoretical 
value based on Eq. (1). The red line indicates a linear fit 
of the data. 

orbit coupling in silicene and hence its topological quan¬ 
tum phase transition. At |e|~0.07, silicene undergoes 
a transition from topological insulator into a metallic 
phase. Within the range of —0.07 <e< -1-0.07, pristine 
silicene remains a QSH insulator with a strain-dependent 
SOC that increases under a compressive strain. The crit¬ 
ical electric field strength is significantly enhanced by the 
screening, nearly 10 times larger than the theoretical pre¬ 
dictions from a tight-binding model. These phenomena 
highlight the interplay between the mechanical strain and 
A—field on the electronic properties of low-buckled hon¬ 
eycomb lattices. 
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